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Abstract 

Small-scale magnetic turbulence observed by the Cluster spacecraft in 
the plasma sheet is investigated by means of a wavelet estimator suitable for 
detecting distinct scaling characteristics even in noisy measurements. The 
spectral estimators used for this purpose are affected by a frequency depen- 
dent bias. The variances of the wavelet coefficients, however, match the 
power-law shaped spectra, which makes the wavelet estimator essentially un- 
biased. These scaling characteristics of the magnetic field data appear to be 
essentially non-steady and intermittent. The scaling properties of bursty bulk 
flow (BBF) and non-BBF associated magnetic fluctuations are analysed with 
the aim of understanding processes of energy transfer between scales. Small- 
scale (~ 0.08 — 0.3 s) magnetic fluctuations having the same scaling index 
a r~j 2.6 as the large-scale (~ 0.7 — 5 s) magnetic fluctuations occur dur- 
ing BBF-associated periods. During non-BBF associated periods the energy 
transfer to small scales is absent, and the large-scale scaling index a ~ 1.7 is 
closer to Kraichnan or Iroshnikov-Kraichnan scalings. The anisotropy char- 
acteristics of magnetic fluctuations show both scale-dependent and scale- 
independent behavior. The former can be partly explained in terms of the 
Goldreich-Sridhar model of MHD turbulence, which leads to the picture 
of Alfvenic turbulence parallel and of eddy turbulence perpendicular to the 
mean magnetic field direction. Nonetheless, other physical mechanisms, 
such as transverse magnetic structures, velocity shears, or boundary effects 
can contribute to the anisotropy characteristics of plasma sheet turbulence. 
The scale-independent features are related to anisotropy characteristics which 
occur during a period of magnetic reconnection and fast tailward flow. 
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1 Introduction 



Magnetic field measurements on board the Geotail and AMPTE/IRM spacecraft 
revealed that the power-law behaviour of the energy density spectra in frequency 
space, P(f) ~ f~ a depends on the considered frequency range, both in the distant 
tail [Hoshino et al, 1994], and in the near-Earth plasma sheet [Bauer et al., 1995]. 
This finding was confirmed by Ohtani et al. [1995] using fractal analysis of the 
fluctuations in the AMPTE/CCE and SCATHA magnetic field measurements. The 
spectral index changes around f c i ~ 0.01 — 0.08 Hz. The observed scaling indices 
are a\ ~ 0.5 — 1.5 for / < / cl and c*2 ~ 1.7 — 2.9 when / > / cl [Milovanov 
et al, 2001; Volwerk et al, 2003, 2004]. The uncertainty of the estimated values 
of a\ and ct2 can be explained by slightly different frequency ranges studied by 
various authors and by the occurrence of spatially inhomogeneous, intermittent 
and transitory mechanisms involved in the generation of the observed scalings. 

Since the characteristic timescale ~ 1/ f c i was found to be a few times the pro- 
ton gyroperiod, ions can play an important role in the excitation of the correspond- 
ing magnetic fluctuations. The high-frequency part of the spectrum can be pro- 
duced by an MHD cascading process. Instability mechanisms or self-organisation 
processes, such as coalescence of magnetic islands or inverse cascades were pro- 
posed for the explanation of low-frequency scaling [Hoshino et al, 1994]. 

For frequencies higher than 0.1 Hz a second break in spectral power density 
was predicted theoretically at a turnover frequency f C 2 [Milovanov et al, 2001]. 
The estimation of the scaling index a 3 for frequencies / > f C 2, however, is en- 
cumbered by noise produced by the magnetometer itself. In this paper we ex- 
amine those high frequency fluctuations for which the spectral indices «2 and 03 
are expected. To this end we use a wavelet method which allows to estimate the 
scaling propeties of magnetic fluctuations over the mentioned frequency ranges. 
Reliable spectral information could not be extracted, hitherto, from the noisy high- 
frequency magnetometer data. Nonetheless, in Section 2 we argue that an essen- 
tially unbiased wavelet method combined with an embedding technique of syn- 
thetic scaling signals allows us to estimate the scaling parameters reliably. 

Very little is known about the origin of the multi-scale turbulence in the plasma 
sheet. Its appearance is obviously dependent on the driving and dissipation mech- 
anisms which characterize an observed fluctuating process [Borovsky and Funsten, 
2003]. Magnetic field and plasma flow fluctuations were observed by a number 
of authors [e.g., Coroniti et al, 1980; Borovsky et al, 1997; Neagu et al, 2002]. 
Plasma shear flows might appear as possible drivers of those fluctuations. The sta- 
tistical analyses of AMPTE/IRM data in the near-Earth magnetotail have already 
shown that there exist rapid plasma flows in the central plasma sheet which are 
preferentially perpendicular to the magnetic field and strongly sunward oriented 
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[Baumjohann et al, 1989, 1990]. Subsequent research indicated that such flows 
are preferentially organized into sporadically occurring groups of bursty bulk flows 
(BBF), typically lasting for 10 min. Despite their short duration, BBFs are the car- 
riers of decisive amounts of mass, momentum and magnetic flux [Angelopoulos et 
al, 1992; Angelopoulos et al, 1993; Schodel et al, 2001]. Therefore BBFs can 
stir the plasma sheet plasma very efficiently. In this paper we place emphasis on 
comparison of multi-scale power and scaling properties of the magnetic fluctua- 
tions during BBF and non-flow or post-flow periods. We will also try to identify 
the turnover frequency f C 2, comparing the estimated values of a 2 and Q3 during 
BBF and non-BBF periods. 

The observed scaling features, together with those expected from MHD consid- 
erations, will help us to classify the magnetic fluctuations observed by the Cluster 
spacecraft in the plasma sheet. Recent approaches to MHD turbulence highlight the 
influence of a mean magnetic field on scale-dependent appearance of anisotropic 
fluctuations [Goldreich and Shridhar, 1995]. The anisotropy increases towards 
small scales. Therefore, the high-resolution Cluster magnetometry provides a way 
for studying the scale-dependent anisotropy of magnetic fluctuations. We discuss 
the observed magnetic anisotropics in terms of recent MHD theories and simula- 
tions. 

2 The wavelet estimator 

We analyse the scaling properties of 67 Hz resolution magnetic field fluctuations 
from the Cluster fluxgate magnetometer [Balogh et al, 2001] using a wavelet 
method proposed by Abry et al, [2000]. Wavelets offer an effective way for get- 
ting an unbiased estimate of a [Veitch and Abry, 1999], whereas power spectral 
estimations are affected by frequency-dependent multiplicative bias [Abry et al, 
1995]. The spectral estimation is biased mainly because it is performed over a 
constant bandwidth which does not match the f~ a structure of the spectrum prop- 
erly. The variance of discrete wavelet coefficients (see later), however, matches 
well the power-law shaped spectrum. Unlike Fourier methods, wavelet methods 
work better for processes exhibiting spectral scaling because wavelets have inbuilt 
rescaling properties. The wavelet estimator allows to identify different kinds of 
scaling processes which are invisible to biased spectral estimators, especially in 
noisy and short data. 
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2.1 Outline of the method 



Abry et al, [2000] proposed a semi-parametric wavelet technique based on a fast 
pyramidal filter bank algorithm for the estimation of scaling parameters Cf and 
a in the relation P(f) ~ Cff~ a , where Cf is a nonzero constant. The algorithm 
consists of several steps. First, a discrete wavelet transform of the data is performed 
over a dyadic grid (scale, time) = (2^,2H) and j, t € N. Then, at each octave 
j = log22 3 ', the variance of the discrete wavelet coefficients d x (j, t) is computed 
through 

^ = ^J2d 2 x (j,t)=2^c f C(a,^) (1) 
n o t =i 

where rij is the number of coefficients at octave j, and C(a, ipo) = J f~ a \ 
rfoif) | 2 df, ipo(f) stands for the power spectrum of the mother wavelet *o- Fi- 
nally, from Equation (1) a and c/ are estimated by constructing a plot of yj = 
log2Hj versus j (logscale diagram, LD) and by using a weighted linear regres- 
sion over the region (j m in, jmax) where yj is a straight line. In this paper we use 
the Daubechies wavelets for which finite data size effects are minimized and the 
number of vanishing moments M can be changed. The latter allows to cancel or 
decrease the effects of linear or polynomial trends and ensures that the wavelet 
details are well defined. 

We note that the wavelet estimator, Equation (1), is also affected by a multi- 
plicative bias (multiplication by C(a, ipo), see later), but C is not frequency depen- 
dent. A similar multiplicative term for spectral estimators is frequency dependent, 
which makes the log-log plot of the power versus frequency not precisely linear. 
Therefore the estimation of a via linear regression in the log-log plot is strongly 
affected by bias [Abry et al, 1995]. Within the wavelet approach Veitch and Abry 
[1999] have found explicit solutions for the variances of log2(t L j) which are then 
used in the weighted regression in the LD. These improvements make the wavelet 
estimator unbiased even for data of small length. 

Different kinds of scaling can be analysed by this technique. The value of a 
and the scaling range (j m ,j n ) are characteristic for the specific scalings present 
[Abry et al, 2000]. For example a € (0, 1) at the largest scales present in the 
data is characteristic for long range dependent processes, while a ~ at large 
scales together with a G (0, 1) at small scales is a signature of noise. Furthermore, 
a G (1,3) aligned at almost all scales available corresponds to a nonstationary, 
self-similar (fractal) process having a fractal dimension D = (5 — a)/ 2 [Abry et 
al, 2000]. 
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2.2 An embedding technique 

When the scaling index is estimated over a high frequency range, magnetometer 
noise can strongly influence the outcome. Since different estimation techniques 
used over slightly different frequency ranges, there is a large scatter in the previ- 
ously observed values of a, which makes it difficult to compare the scaling prop- 
erties in turbulence. Furthermore, the estimated values are presumably affected by 
the occurrence of transitory processes in the plasma sheet. We therefore perform 
the estimation of Cf and a within sliding overlapping windows of width W with 
a time shift S < < W to analyse the data over different frequency ranges with the 
same method and simultaneously take into account their temporal evolution. 

In order to examine how the signal-to-noise ratio (SNR, the relative power of 
the signal and the noise in LD) influences the estimation of a within the window 
W in small-scale fluctuations, we generate a self-similar fractal signal with known 
second order scaling properties using wavelets and embed it in a time series of quiet 
time magnetic field measurements from the Cluster spacecraft. As a synthetic frac- 
tal signal we use fractional Brownian motion (fBm) generated by a fast filter bank 
based pyramidal algorithm [Abry and Sellan, 1996]. Fractional Brownian motion 
has a Gaussian distribution of increments by definition. At the same time, mag- 
netic fluctuations exhibit two-point Gaussian statistics only for large separations. 
Therefore, we consider the fBm process as an approximation which can roughly 
model the second order properties of the magnetic fluctuations. 

We suppose that the magnetic fluctuations observed in the plasma sheet during 
quiet times are close to the magnetometer noise level. Figure la shows the mag- 
netic field Bx component (67 Hz resolution) observed by Cluster 3 on Septem- 
ber 7, 2001. If not specified the GSM coordinate system will be used through- 
out the paper. Two short fractal signals A and B, with the same scaling index 
ola = oib = 1.6 but different amplitude are added to Bx- In Figure la Bx is 
depicted from 1840 to 1917 UT. First, signal A is added into Bx from ~1857 to 
~1901 UT. The average amplitude of fluctuations of A is much larger than of Bx- 
In Figure lb a dashed horizontal line show the theoretical value of ola = «b = 1-6. 
The estimated a a reaches this theoretical value in a wider interval in Figure lb as 
was the length of signal A in Figure la. Despite the non-precise timing, the es- 
timated value of a a is correct. Next, signal B is embedded into Bx- Signal B 
and Bx have comparable average amplitude. The estimated as in Figure lb does 
not achieve the theoretical horizontal line, therefore ag is underestimated. We can 
conclude that, whenever the amplitude of the signal is significantly larger than the 
amplitude of the noise (case A), the correct value of a is recovered, otherwise a is 
underestimated (case B). 

Let us explain now the reason of non-precise timing. In Figure lb as are esti- 
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mated within a sliding window W ~ 480 s over time scales (octaves) j\ = 2 and 
j 2 = 4 (~ 0.08 — 0.33 s and ~ 3 — 12 Hz; these are approximate Fourier peri- 
ods and frequencies that correspond to the main oscillations within the Daubechies 
wavelets). The embedded fractal signals A and B produce a rapid departure of the 
estimated values a from the noise level a ~ 0.3, as soon as the analysing window 
includes some part of a fractal signal (A or B in Figure la). This effect is visi- 
ble at ~1853 UT and before 1906 UT, when two analysing windows depicted by 
horizontal lines in Figure lb are centered in those time positions. The dashed, ver- 
tically extended right and left ends of the windows match perfectly the beginnings 
and the ends of the fractal signals A and B. When the timing is important, the effect 
of window length has to be accounted. 

Another effect which influences the estimation of the scaling index a. If the 
duration of a scaling process T is less than half the window length, W/2, a is also 
underestimated. In Figure la the length of the fractal signals A and B equals W/2. 
The proper sliding window length for unbiased estimation of a can be found by 
trial. For decreasing W the estimated a increases, then saturates. Obviously, as 
W decreases, the fluctuations of the estimations increase. In general, W has to 
be chosen such that the opposing requirements for robustness of scaling parameter 
estimations (wide window needed) and for time localization of short events (narrow 
window needed) were matched. For the events analysed in this paper W = 61 s 
was found to be an optimal choice. 

Finally, let us consider possible kinds of non-physical mechanisms which can 
lead to spurious scalings. Time series analysis methods may lead to spurious results 
when applied to short data series containing irregularities, e.g. jumps, spikes [Kat- 
sev and L'Heureux, 2003] or random pulse trains [Watkins et al, 2001]. Spikes, 
i.e., single values with large deviation are removed during quality check of the data. 
Possible effects of other irregularities will be discussed jointly with the estimation 
of the scaling parameters in subsequent sections. 

2.3 Noise correction 

The above explained embedding technique allows us to introduce corrections, when- 
ever the amplitude of a fractal signal approaches the noise level. It can happen 
when high frequency (small-scale) scalings are examined. A few more steps are 
needed to explain the method in detail. To this end we show the typical scaling 
properties of a fractal signal (s), noise (n) and of a compound additive process 
(signal+noise, sn) in LD. For what follows the subscripts s, n, sn will correspond 
to the specified processes. We also demonstrate that a single correction curve (see 
later) belongs to every embedded model fractal signal with a given a s . 

Equation 1 describes the scaling properties of a signal through the scaling of 
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the corresponding wavelet coefficients. Taking a base 2 logarithm of both sides of 
Equation 1, we obtain an equation of line in LD. Therefore any scaling signal (a 
straight line in LD) can be characterized by a pair of values (c/(j), a(j m in,jmax)), 
where Cf(j) is the power of intercept at octave j. Here we suppose that, for a fixed 
mother wavelet, the last term, C(a,ipo), on the right hand side of Equation 1, is 
dynamically unimportant. The change of i/jq, however, can influence the estimation 
of scaling parameters. We use the Daubechies family of wavelets with changing 
number of vanishing moments M. It does slightly influence our estimations and 
the corresponding uncertainities represent a source of error. 

Figure 2 shows an exemple of LD. The magnetometer noise is depicted by a 
single curve, since we suppose, it is time-independent (see later). It exhibits signa- 
tures typical for noise: a n ~ at large scales (j > 7) together with a n ~ 0.3 € 
(0, 1) at small scales (j < 6) [Abry et al. 2000]. For the sake of clarity the straight 
lines corresponding to the scalings of both signal and compound process are de- 
picted only at the left hand side of Figure 2 (j < 5). The amplitude of the model 
fractal signal was changed from values larger to values smaller than the amplitude 
of the noise. The resulting parallel, shifted straight lines are representations of self- 
similar signals in LD. The compound process exhibits similar characteristics as the 
signal itself, while SNR is large. If SNR becomes < 1, the compound process 
approaches the straight line of noise at the scales j± and j2- 

Time-independent noise has averages ((c/ n (j)}, (a n (ji, ji))) ~ const. Our 
goal is to estimate the time evolution of the scaling exponent of the signal a s (ji ,]2,t) 
for known (c /n , a n ) and (c' fsn (j, t), a sn {ji, j 2 , t)). The normalization c' fsn (j, t) = 
c fsn(j,t)/(cf n (j)} ensures that for weak signals, Cf sn ~ Cf n (noise dominates), 
and Cf sn ~ 1. Actually, Cj sn (t) represents a time dependent SNR. 

Now we can use the introduced notation and the embedding technique to study 
the dependence of a sn on c^ sn within a window W. Figure 3 shows the relation- 
ship a sn (c'j- sn ) (correction curves), which was obtained experimentally by embed- 
ding two model fractal signals with a mo d = 1-5 and 2, gradually decreasing their 
amplitudes (decreasing SNR), into quiet time Bx measured by s/c 4 on Septem- 
ber 7, 2001. For both model fractals the correction curves were computed at the 
octaves j = and j = 4, separately. Independent of j, a sn approaches asymptot- 
ically a s = a mo d for large enough Cj sn . Large Cj sn means that SNR is large and 
the compound signal exhibits the same scaling as the fractal model signal (left-top 
straight lines in Figure 2). If Cj sn ~ 1, the noise dominates, and a sn ~ a n (scaling 
of the compound process between yj = —5 and —10 in Figure 2). Between the 
two limits (cj sn — > oo and Cj sn ~ 1), however, the course of the curves for j = 
and j = 4 in Figure 3 is rather different. Because of the complex behaviour of the 
compound process at j = near the noise level (Figure 2), a sn cannot be deter- 



7 



mined from the a sn (c'j sn ) relationship unequivocally (Figure 3). Since the relation 
is unambiguous at the scale j = 4, we will consider the a sn {c^ sn ) curves at that 
octave throughout the paper. 

The accuracy of the wavelet method estimation of a can be controlled by 
the power of intercept of the compound process c^ sn . Figure 3 shows that for 

c 'fsn we fi 11 ^ a sn < ce mo d, therefore a s is underestimated. The a sn (c'f sn ) 

relationship allows us to introduce corrections and estimate the effective asymp- 
totic value of a s . 

We note that Cf sn (t) and a sn (ji,j2,t) are computed directly from magnetic 
field time series by using a sliding window W. The key assumption of the method 
is time-independence of the scaling parameters for noise. Obviously, this condition 
is not strictly fulfilled. Therefore, it is important to examine the time-independence 
of (cf n ) and (a n ). The parameters Cj sn (i) and a sn (ji,j2, t) computed from mag- 
netic field Bx component (67 Hz resolution) observed by the Cluster spacecraft 
s/c 1 and 3 during the quiet interval 1700-2100 UT on September 7, 2001, are de- 
picted in Figure 4. Here W = 68 s and S = 4 s. During the chosen interval the 
spacecraft are traversing the neutral sheet from the northern {Bx ~ 18 nT) to the 
southern hemisphere {Bx ~ —15 nT) (Figure 4a) and the plasma flow velocity 
does not exceed 50 km/s (not shown). Independent of the location of the Cluster 
spacecraft and despite the occurence of low frequency wave motion after 1830 UT, 
both a sn (j\,j2,t) and Cj sn (i) exhibit stationary behaviour, fluctuating around the 
mean values (Figure 4 b,c). The estimated a sn can be related to the magnetome- 
ter noise a n because a sn ~ a n <G (0, 1) and Cj sn ~ Cf n ~ 1 (see also Figure 
3). Noteworthy is that a n and Cf n are slightly different for each s/c (Figure 4b,c) 
and also vary with the magnetic field components. The examination of several 
quiet time magnetic field time series and time series obtained from ground-based 
calibration measurements shows that in average (including the four s/c and each 
magnetic component) {a n ) = 0.5 ± 0.25 and (cf n ) = 0.004 ± 0.003. Similar 
values were obtained using 128 Hz quiet time magnetic field measurements from 
the Equator-S spacecraft (not shown). For discrimination between a s (which is 
expected to be > 1) and a n , this accuracy seems to be sufficient. It is expedient 
in concrete case studies, however, to compute the correction curves (Figure 3) for 
each spacecraft and component separately. 

Examining a sn and Cj sn estimated for Bx from s/c 2 on September 7, 2001, 
we found short, but significant deviations from (a n ) and (cf n ), respectively (not 
shown). The deviations occurred due to a 5 min long data gap on s/c 2, which 
led to spurious scalings because of the step function-like irregularity in Bx, in ac- 
cordance with the findings of Katsev and L'Heureux [2003]. This type of spurious 
estimations can easily be avoided by interpolation. Another type of jump-related 
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spurious estimation might occur near the plasma sheet boundary layer (PSBL) or 
owing to its motion. The PSBL cannot be identified from the magnetic field mea- 
surements alone [Eastman et al, 1984]. Therefore, non-physical scalings due to 
magnetic jumps near PSBL are not expected. 

In summary, using the wavelet approach, we can estimate directly from a time 
series the time evolution of the parameters of a compound process (cj sra (t) , a sn (t) ) . 
The time independent parameters of noise ( (c/ n ) , (a n ) ) can be obtained from quiet 
time measurements. When Cj sn (t) is small, the effective asymptotic value of a s (t) 
for a self-similar signal can be obtained from a correction curve taking a limit 
Cf sn (t) — ► oo. Spurious estimations during quiescent periods can be avoided by 
careful quality checks of the magnetic data. 

3 Magnetic fluctuations in the plasma sheet 

In order to introduce the wavelet method and the small-scale corrections to the esti- 
mated a, quiet time magnetic measurements together with artificial fractal signals 
were used in the previous section. Since our goal is to investigate BBF and non- 
BBF related magnetic turbulence, here we choose a ~ 4h interval during which the 
Cluster spacecraft observed BBF and non-BBF intervals successively. 

3.1 Event overview and small-scale corrections 

In what follows, we analyse burst mode (67 Hz) magnetic data from the Cluster 
fluxgate magnetometer (FGM) [Balogh et al, 2001] during the interval 0110-0500 
UT on August 27, 2001, when the Cluster spacecraft were near the Zqsm = 
plane, in the postmidnight magnetotail (Xqsm ~ — 19-Re). The magnetic data are 
compared with the spin-resolution (4 s) velocity data from the Cluster ion spec- 
trometry (CIS/CODIF) experiment [Reme et al, 2001]. Figure 5a shows the Bx 
component from s/c 1,3. The variation of the scaling index for two different ranges 
of scales is depicted in Figure 5b. In addition to (ji, j'2) = (2,4), corresponding 
time scale ~ 0.08 — 0.33 s), we also consider the octaves: (j'3,,7'4) = (5, 8), corre- 
sponding time scale ~ 0.7 — 5 s . Hereinafter we will refer to those scales as small 
scales (ji, j'2) and large scales (j'3, J4), respectively. Obviously, j\ is limited by 
the resolution and j'4 by the chosen window W. The average noise levels (a n ) for 
s/c 1,3 and Bx magnetic components are also indicated by black and gray dashed 
lines in Figure 5b. The large-scale scaling exponents are anticorrelated with the 
magnitude of Bx simply because of the preferential occurrence of perpendicular 
flows closer to the neutral sheet (Bx — > nT). Figure 5c shows the variation of 
the normalized power of the intercepts c fsn computed for the small scales (ji,^)- 
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Its time evolution exhibits periods of c'j sn ~ 1 intermittently interrupted by bursts 
of activity when c'j sn > 1. The former corresponds to the noise level, indicating 
that there is no significant energy flow from large scales into small scales. The 
increased values of Cj sn are associated with high frequency (small-scale) fluctua- 
tions of Bx (Figure 5a), increased a sn (relative to (a n ), Figure 5b) and are well 
correlated with the appearance of bursty perpendicular flows (Figure 5d). The large 
scales (j3, j'4) are not corrupted significantly by magnetometer noise. 

The effect of spurious estimations due to spikes, random pulses, or other types 
of nonstationary signals can be addressed for BBF-associated periods, too. The 
use of Daubechies analysing wavelets with vanishing moments ensures that part of 
the nonstationarities represented by polynomial trends are cancelled. Spurious es- 
timates of the scaling index might appear because of discontinuities which separate 
earthward moving plasma structures from the background plasma ahead of them 
[Sergeev et al, 1996]. In fact, before the appearance of the rapid perpendicular 
flow, e.g., at 0131 UT (Figure 5d), there is a large jump in the large scale scaling 
index from a sn (js, j'4) ~ 1 to 2.5. This jump and all similar jumps in our data can 
partly be explained by the presence of a discontinuity in front of the plasma flow 
[Nakamura et al, 2002]. On the other side, this short effect cannot be separated 
from the non-physical effect of the sliding window, explained above (Figure 1). 

After 0130 UT, a sn (jz, j±) fluctuates around ~ 2.5 for more than 30 min (Fig- 
ure 5b). We believe, that all estimations exhibiting small fluctuations around a 
mean value for equal or longer time than the sliding window length (~ 1 min), 
belong to a BBF-associated scaling process in our case. Therefore these estima- 
tions are not spurious. The same is true for non-BBF associated processes. Also, 
a sn is largely coherent between s/cl and 3. Short living pulses are supposed to 
be random, however. Moreover, some results of multifractal approach [Voros et 
al., 2003] and of wavelet based structure function approach [Voros et al, 2004] 
also support the idea that BBF-associated magnetic fluctuations correspond to a 
physical scaling process. 

During the BBF-associated periods one finds that 1 < Cj sn < 100, therefore 
a s is underestimated (see Figure 3). Since self-similarity is defined for 1 < a < 3, 
we used 4 model fractal time series with a mo d = 1.5,2,2.5,3 and quiet time 
Bx magnetic data from s/c 1,3 to compute the a sn (c^ sn ) correction curves at the 
octave j=4 (Figure 6). These curves can be used to introduce corrections and to 
estimate the effective value of a s . Like in Figure 3 there are two limiting cases. 
The condition c^ sn — > 00 leads to a sn = a s = a mo d, with slower approach to 
a model value for a sn > 2. The limit c^ sn — > 1 leads to a sn = (a n ). The 
error bars depicted at the left-top corner of Figure 6 arise from the uncertainty 
of the estimation of a sn in a logscale diagram using Daubechies wavelets with 
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varying number of vanishing moments (M = 3 — 7) and from the uncertainty of the 
magnetometer noise level estimation within the chosen window W. The depicted 
error bars represent only the uncertainity of the determination of the correction 
curves, and not the uncertainity of the estimations of a s (j\, ji). The estimation 
error of the small scale a changes from ±0.2 to ±0.5 as Cj sn — > 1. The errors can 
be estimated by adding more correction curves to Figure 6 (not shown). It is more 
difficult to distinguish an a as c^ sn decreases. We present two examples how to 
read the measured a sn (ji, j-i) using both Figures 5 and 6. The two dashed lines 
A and B in Figure 6 correspond to two pairs of values (cj sn , a sn ) in Figure 5b,c 
before 0330 UT. The gray lines and the letter A belong to s/c 3, while the black lines 
and the letter B belong to s/c 1. In both examples the corrected effective value of 
small-scale a s (ji, j'2) is close to the correction curves which asymptotically reach 
the horizonthal line a mo d = 3 in Figure 6. After correction, all BBF associated 
peaks of a s (j\ , j'2) > (a n ) and Cj sn > 1 show a similar feature (not shown), that is 
in average, a s (ji,j2) = 2.5±0.5 ~ as(j3, J4) ~ 2.5. Therefore, after introducing 
corrections, the BBF-associated small-scale and large-scale scaling indices cannot 
be discriminated. We note that the peaks at ~ 0240 UT from s/c 1 (Figure 5b,c) 
appear due to the data gaps. 

3.2 Statistical evaluation of BBF and non-BBF associated multi-scale 
magnetic fluctuations 

During non-BBF periods a s (js, j'4) drops to the noise level when s/c 1,3 are in the 
lobe (Bx ~ 25nT) or fluctuates between 1 and 2.5 when s/c 1,3 are in the plasma 
sheet (Figure 5a,b). Among others, the time evolution of the scaling indices can be 
influenced by the location of the Cluster spacecraft, the duration and the magnitude 
of perpendicular flows and by the motion of the plasma sheet (e.g., thickening or 
thinning). Despite the complexity, we can separate BBF and non-BBF associated 
magnetic fluctuations, based on two different conditions. 

Condition 1 is based only on small-scale scaling parameters having for BBF- 
associated fluctuations on both s/c 1,3, a sn (ji,j2) > (a n ) ±0.3 and c^ sn > 2 and 
non-BBF-associated fluctuations a sn (ji, ji) < (a n ) ± 0.1 and Cj sn < 1.1. It sim- 
ply means that we relate the observed small-scale magnetic fluctuation-dissipation 
activity to the large-scale magnetic fluctuations. We suppose that on the basis of 
scaling characteristics and power of small-scale magnetic fluctuations we can iden- 
tify BBF and non-BBF associated intervals without actual examination of velocity 
data. For testing this supposition we use only the parameters computed from the 
magnetic field Bx component measured by s/c 1,3 (Figure 5a-c). So we select 
only those values of a s (js, j'4) from the whole period (Figure 5b) which meet the 
above conditions for a sn (ji, J2) and Cj sn . The conditional statistics performed on 
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the large-scale scaling index gives for both s/c 1,3 a s ,BBF(j3, ji) = 2.57 ± 0.05 
for BBF-associated and o. s ,n<m-BBF{h-, 3a) = 1-6 ± 0.4 for non-BBF associated 
case. 

Now we consider also the velocity data. Condition 2 is based on threshold- 
ing of perpendicular flows, taking V± jXy = \JV± X + V± y > 300 km/s for rapid 
flows and V± yXy < 100 km/s for non-BBF associated remnant flows. Since the 
velocity data is available with a time resolution of ~ 4s and the time step between 
neighbouring sliding windows is S = 4s, after adjusting times, a s (js, ji) can be 
compared to V± jXy . It gives a SjB BF{jz,k) = 2.6±0.1 and a Stnon - BBF (j 3 , j 4 ) = 
1.7 ± 0.4. All the estimates represent mean values with standard deviations. 

Both conditions give the same statistical results on large-scale scaling index. 
From the conditions and the corresponding values of large-scale indices follows 
that, relative to non-BBF associated intervals, BBF-associated intervals exhibit 
both large-scale and small-scale spectral steepening and an unambiguous increase 
of small-scale power of magnetic fluctuations. Therefore we propose that small- 
scale magnetic fluctuations are energized by BBF's. Then the spectral transfer of 
energy to small scales is transient, which makes the position of the spectral break at 
the expected turnover frequency f C 2 uncertain. Our results suggest that whenever 
BBF-associated magnetic fluctuations occur, the break at ~ f C 2 disappears when 
the small-scale corrections are introduced. Otherwise the energy flow to small 
scales ends at f C 2, and in this case, at frequencies / > f C 2, the measured scaling 
parameters belong to the noise. The examination of the LD's at various ranges of 
scales shows that f C 2 may correspond to the approximate Fourier frequencies of 
1 — 3 Hz (between octaves j2 and j'3). 

The same analysis of the By and Bz magnetic field components shows that 
the observed scaling features remain similar to the Bx case. The small-scale 
By and Bz fluctuations are also associated mainly with V± :Xy . After introduc- 
ing corrections as above, a s (ji, j'2) ~ ot s (jz,j±) ~ 2.6 for BBF and ~ 1.7 for 
remnant flows, the errors in the same range as above, in conditional statistics of 
Bx- A dissimilarity appears in Bz small-scale fluctuations (Figure 7). The small- 
scale power Cj sn of Bz fluctuations (Figure 7c) during BBF-associated intervals 
is several times larger than that of Bx (Figure 5c) or By (not shown). Usually 
Cf sn (Bz) » Cj sn (By) > Cf sn (Bx) over time scales less than 1 s. This is a 
signature of anisotropy and we will examine it deeper. 

3.3 Anisotropy of magnetic fluctuations 

The small-scale powers of the magnetic field components from s/c 1,3, computed 
at j = 4 (~ 0.33 s), are compared in log-log scatterplots in Figure 8. The 
comparison of Bx and By components shows a balanced scatter around the line 
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Cf sn {Bx) = Cf sn {By) (Figure 8a), though the increased population of points un- 
der that line between 2 < c^ sn (By) < 10, indicates a weak anisotropy. The scat- 
terplot of c'f sn (Bx) and c'j- sn (Bz) shows strong anisotropy (Figure 8b) for values 

c 'fsn > 1-5- We also calculated the scatterplots at scale j = 6 (~ 1.3 s) (Figures 9 
a, b), which, after comparison with Figure 8b, indicate scale dependent anisotropy. 
The anisotropy increases towards small scales. Since Cj sn (Bz) S> Cf sn (Bx) and 

Cf sn {By) > c'f sn (Bx) in Figure 8, the small-scale magnetic fluctuations are gen- 
erated due to multi-scale transfer of energy preferentially in the GSM Z or Z — Y 
direction. The powers at the scale j = 6 show more isotropy, even if the scatterplot 
of Cf sn (Bx) versus Cf sn (By) (Figure 9a) indicates that, with increasing scales the 
field-aligned power of fluctuations can also increase. 

The points depicted by '+' signs in Figures 8b and 9b belong to the time period 
between 0400 and 0410 UT. During this interval a magnetic reconnection event 
in the central plasma sheet at 0401 UT was identified by Baker et al. [2002]. 
The corresponding '+' points in Figures 8b and 9b, therefore, are associated with 
the magnetic reconnection or/and the tailward flow, leading to a strong anisotropy 
which does not exhibit the same increase towards the small scales as it is observed 
for the other points. Therefore the underlying mechanism which drives these fluc- 
tuations should be different from the previous scale dependent one. Since only this 
reconnection-related short interval shows different anisotropy characteristics, we 
cannot say more about the nature of the underlying physical processes. Answering 
the question about reconnection or tailward flow associated anisotropy character- 
istics certainly requires a thorough statistical analysis. We can speculate, however, 
that the observed scale-independent anisotropy features might be related to the spe- 
cific current systems or magnetic field topology near the reconnection site [Nagai 
et al, 2003; Runov et al, 2003]. 

The small-scale a sn (ji, j-i) for individual components are compared in scatter- 
plots in Figure 10a,b. Again, a sn (ji,j2,Bz) seems to be different from a 8n (ji,j2, Bx) 
and a sn (ji,j2, By)- This anisotropy can weaken, however, because non-corrected 
values are depicted and the corrections are larger for Bx , By than for Bz (Figures 
5c and 7c). 

The anisotropy can also be considered in the coordinate system of the local 
mean field. We computed the mean total magnetic field vector through low-pass 
filtering the GSM magnetic field components and using standard vector operations 
transformed the burst mode magnetic data into mean-field perpendicular (B±) and 
mean-field aligned (B\\) components. Figures 11a, b show B± and B\\, respec- 
tively, computed for averaged (B) at the time-scale 1.3 s. Figures 11 c, d show the 
perpendicular Cf sn (B±) and the parallel Cj sn {B\\) powers, computed at the scales 
j = 3 (~ 0.33 s) and j = 6 (~ 1.3 s). As above for the GSM X, Y, Z components 
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(Figures 8-10), the anisotropy is scale dependent and the small-scale power pref- 
erentially develops in directions perpendicular to the mean magnetic field. Again, 
the short interval after 0400 UT shows a different kind of reconnection or tailward 
flow related behavior. The large-scale (j = 6) perpendicular power remains larger 
than the parallel power only during this interval. The similarity of the statistical 
results obtained in GSM and mean field coordinate systems can be explained as 
follows: for open, stretched magnetic field lines B\\ resembles Bx and B± is sim- 
ilar to yj B\ + By (where By and Bz are the GSM components of the magnetic 
field). Dissimilarities can appear close to the neutral sheet when Bx — > 0. 

3.4 Two-point conditional statistics 

Let us investigate now the effect of the transitory character of the fluctuations or 
their drivers. We compare magnetic fluctuations observed at the same time on two 
spatially separated Cluster spacecraft. Figure 12 compares the parallel a s \\(jz,ji) 
and the perpendicular a s ±(j3, Ja) scaling indices from s/c 1 and 3. The scatterplot 
in (Figure 12a) involves all the magnetic data from s/c 1,3 and forms an irregu- 
lar cloud of points. We applied two-point conditional statistics to find BBF- and 
non-BBF-associated points in that cloud. We used the same small-scale condi- 
tion 1 combined with the velocity condition 2 as above (both should give the same 
results), and required the fulfilment jointly at s/c 1,3 at the same time. In conse- 
quence of the joint conditional requirement the amount of data points in (Figures 
12 b,c) is significantly reduced. It can be explained partially by a relatively high 
level of fluctuations of the considered parameters. For example, the perpendicular 
velocity wildly fluctuates during the rapid flows (e.g. Figure 5d) which makes the 
comparison with the small-scale parameters more difficult. The omission of the 
velocity condition, however, does not add many more points to the Figures 12 b,c. 
Therefore, we suppose that fewer joint data points are found by two-point condi- 
tional statistics because the two spacecraft are often in physically different regions, 
exhibiting different scaling features. It indicates that the magnetic fluctuations and 
their sources are transient and limited in space. 

The statistical evaluation of the BBF-associated data points in Figure 12b gives 
averages a s ±(j3, j±) = 2.4 ± 0.3 and ck s || (j3, Ja) = 2.7 ± 0.3 which is consistent 
with the previous estimations using the components Bx,y,z (dashed lines). The 
'+' signs in Figure 12b represent measurements during the reconnection event, 
discussed above. 

Figure 12c shows the non-BBF-associated data points. Here the '+' signs rep- 
resent measurements near the plasma sheet boundary layer or in the lobe (Bx > 20 
nT) while the remaining points are inside of the plasma sheet (Bx < 15 nT). The 
statistical evaluation of the latter gives a s ±(js, j'4) = 1.7 ± 0.3 and a a ii(j3, j'4) = 
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1.7 ± 0.4, which is again consistent with the previous estimations using the com- 
ponents Bx,y,z (dashed lines). 

Let us summarize the main results at the end of this section. We used a wavelet 
method to study the small-scale properties of magnetic fluctuations. We found that 
the energy transfer to the smallest available time scales is associated with BBFs. 
The scaling properties of those fluctuations are recoverable from burst mode mag- 
netic fluctuations by using an embedding technique. Small-scale fluctuations show 
an increased anisotropy relative to the mean magnetic field. The anisotropy is ob- 
served mainly in the normalized power of the magnetic fluctuations, while it is less 
evident in the spectral scaling index. The transient character of the physical pro- 
cesses or driving mechanisms makes a robust estimation of the basic characteristics 
of magnetic turbulence difficult in the plasma sheet. 

4 Discussion 

A complete characterization of MHD turbulence requires not only a statistical de- 
scription of the velocity field, but also a statistical study of other measurable quan- 
tities, e.g. magnetic field, density or temperature fluctuations. In this paper we 
analysed a 4 hour long period of high quality burst mode (67 Hz) magnetic field 
data, which allowed an immersion to the scale or frequency ranges which are not 
available in velocity measurements. We considered time scales from ~ 0.08 to 5 
s, which, assuming 500 km/s velocities during rapid flows, correspond roughly to 
spatial scales of 40 — 2500 km in the flow direction. The average velocities are 
smaller than 500 km/s, so the largest scale can be significantly smaller than 2500 
km. Also, both the velocities and the estimated spatial scales are smaller in other 
than the maximum velocity direction or during non-BBF intervals. Since the av- 
erage separation of the Cluster spacecraft was ~ 1500 km during 2001, we used 
principally a single-point technique to analyse scaling characteristics of magnetic 
fluctuations from s/c 1,3. The usual multi-point statistical approach, commonly 
available in laboratory measurements, is not available in the plasma sheet at those 
scales and the conversion of the scaling characteristics of magnetic fluctuations in 
time into scaling laws in wavenumber space is possible only through additional hy- 
potheses. The Taylor hypothesis assumes that the conversion is possible if the spa- 
tial fluctuations on a given scale pass over the spacecraft faster than they typically 
evolve in time. In the plasma sheet this can be the case during fast BBF's [Hor- 
bury, 2000] or during strong global bulk flows which may occur after a dynamical 
reconfiguration of the magnetotail [Borovsky and Funsten, 2003]. Otherwise, the 
random velocities of the fluctuations advect spatial structures over the spacecraft. 
A more realistic approach is the so-called random sweeping model (RSM), which 
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takes the basically random velocities of the flow into account and allows to obtain 
statistical information on spatial structures. 

Borovsky et al. [1997] have used RSM and found that the spectral index, de- 
scribing the scaling properties of magnetic fluctuations in the plasma sheet at the 
frequency range of 1/4 — 1/150 is the same in both frequency or wavenumber 
representations, that is, k~ a ~ f~ a . The 'large-scale' (~ 1/0.3 — 1/5 s" 1 ) used in 
this paper partially overlaps with the scales of Borovsky et al. [1997] paper. Note 
that, RSM is supposed to work when the large scale velocity field drives the small- 
scale fluctuations. This seems to be valid in our case due to the high correlation 
between BBFs and small-scale magnetic fluctuations. Therefore we suppose, that 
at least at the scales (j'3, j'4) the equality k~ a ~ f~ a holds also in our case and the 
estimated values of a can be compared with the power-law dependence of power 
spectra in MHD turbulence. We have to underline, however, that RSM is only a 
model and its validity should be proven comparing real spatial and time structures. 
In absence of such comparisons the interpretation of scaling indices in terms of 
MHD turbulence remains disputable. At the same time, we consider the observed 
anisotropy effects to be more reliable. Let us start the comparison of observations 
and theory, therefore, with the evaluation of the influence of a local mean magnetic 
field on anisotropy of fluctuations. 

The clear signs of the interaction of rapid flows with the magnetic field in 
the plasma sheet are the BBF-associated dipolarization of field lines (increase of 
Bz) and the breaking and diversion of the flows around the dipolar inner magne- 
tosphere. The presence of a mean magnetic field in the plasma sheet cannot be 
eliminated by coordinate transformations. This is in contrast with hydrodynamic 
turbulence, where the anisotropy usually ceases in the coordinate system of the av- 
erage flow [Biskamp, 2003]. The classical Kolmogorov theory for hydrodynamic 
and the Iroshnikov-Kraichnan theory for MHD turbulence assume isotropy of the 
inertial-range energy cascade in Fourier space, exhibiting scalings E{k) ~ /c~ 5 / 3 
and ~ k~ 3 / 2 , respectively (E(k) is the amount of energy between wavenumbers k 
and k + dk divided by dk). 

Shebalin et al. [1983] studied incompressible MHD anisotropics arising in 
wave vector space in the presence of a mean magnetic field. They studied the in- 
teraction of opposite-travelling wave packets and found that, in wave vector space, 
those interactions produce modes with wavevectors preferentially perpendicular to 
the mean magnetic field. Goldreich and Sridhar [1995] proposed a balance condi- 
tion between parallel and perpendicular modes (parallel propagating Alfven waves 
and perpendicular eddy motions). On this basis it was shown that a scale dependent 
anisotropy appears, namely fey ~ k 2 ^ 3 (or /y ~ I 2 / 3 , where I ~ l/k is the char- 
acteristic scale of eddies in parallel and perpendicular directions), which means 
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that the anisotropy is increasing with decreasing scale. Moreover, it was shown 
that the perpendicular fluctuations exhibit Kolmogorov scaling E(k) ~ k~ 5 / 3 . In 
other words, the mean magnetic field does not influence the perpendicular hydrody- 
namic cascade formed by eddies. The coupled parallel fluctuations are essentially 
produced by waves. In the case of compressible MHD the Alfven modes follow 
the anisotropic Goldreich-Sridhar scaling while fast modes exhibit isotropy [Cho 
et al, 2003]. 

The scale dependent anisotropy proposed by Goldreich and Sridhar [1995] 
would represent a way the relative importance of Alfven waves and eddy motions 
might be recognized. Matthaeus et al. [1998] tested this assumption using Fourier 
spectral methods for MHD equations in a periodic cube. They have found, that 
for incompressible, weakly compressible and driven MHD turbulence, anisotropy 
scales linearly with the ratio of fluctuating to total magnetic field strength. Their re- 
sults therefore seem to be inconsistent with the Goldreich-Sridhar model. Cho and 
Vishniac [2000] argued, however, that the Fourier technique used by Matthaeus 
et el. [1998] smooths out the true scaling relation for anisotropy. It happens be- 
cause the anisotropic eddies or wave packets are elongated along magnetic field 
lines. When they follow the large-scale variation of the magnetic field, they point 
in different directions. Consequently, the information about the eddy shapes is lost 
when a global transformation such as the Fourier transformation is applied. Cho 
and Vishniac [2000] showed that, when the eddy shapes relative to the local mag- 
netic field are analysed in real space, the Goldreich-Sridhar scaling is recovered. 

For what follows, it is important that not only the velocity statistics exhibits the 
Goldreich-Sridhar anisotropy, but also the magnetic field statistics [Cho and Lazar- 
ian, 2004]. If we suppose similarity between frequency and wavenumber scalings, 
then our results suggest that the observed magnetic turbulence in the plasma sheet 
also resembles the anisotropy of Goldreich-Sridhar model. The small-scale mag- 
netic turbulence is stronger in perpendicular directions. It is mainly visible when 
the relative powers of fluctuations are compared (Figures 8 and 11). Then, ex- 
cept for the short interval when magnetic reconnection and tailward flow occurs, 
the anisotropy is scale dependent (Figures 8,9). The applied wavelet technique 
ensures that the scale-dependent anisotropy is not smoothed out. Namely, over 
each considered scale range the local trend (local mean field) is cancelled by the 
analysing wavelet having a proper number of vanishing moments. The observed 
anisotropy seems to indicate that the small-scale eddy turbulence is stronger than 
the Alfvenic turbulence. This is in agreement with the results of Borovsky and 
Funsten [2003], who claimed that turbulence of eddies rather than the turbulence 
of Alfven waves prevails in the plasma sheet. Their observation is based on three 
different arguments: (a.) Alfven wave turbulence produces scale sizes that are 
larger than the vertical size of the plasma sheet, whereas eddy turbulence does not; 
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(b.) the autocorrelation length of the turbulence is typical for eddy turbulence; (c.) 
the low Alfven ratio (ratio of the kinetic energy of the flow to the energy den- 
sity of magnetic fluctuations) indicates 2D eddy structures. All these results might 
yield a picture of the Goldreich-Sridhar model and quasi-2D MHD turbulence. The 
quasi-2D MHD picture is supported by numerical studies of high Reynolds number 
turbulence which predict that the spatial distribution of the turbulent fields is more 
intermittent in 2D than in 3D MHD ([Biskamp, 2003], and references therein). 
Since BBF-associated magnetic fluctuations are stretched down to the small scales 
(jij 32) (Figures 5,7), where the anisotropy becomes stronger and turbulence might 
be essentially 2D, an increase of the small scale intermittence can be expected, 
too. Indeed, Voros et al. [2003] have shown that small scale magnetic field in- 
termittence increases during BBF-associated intervals. Here, some more words of 
caution have to be added. First of all, the importance of multi-scale dynamics in 
turbulence has to be underlined. Since the Goldreich-Sridhar anisotropy is scale 
dependent, at sufficiently large scales this kind of anisotropy can disappear, and 
even Alfvenic turbulence can become stronger. Figures 9a and lid show that the 
large scale power of parallel (Bx or B» ) fluctuations can be equal with or slightly 
overpower the perpendicular components (By, Bz or B±). The development of 
anisotropy and 2D structures can also be dependent on the actual Reynolds num- 
bers, helicity, Alfven ratio [Oughton et al, 1994] or be influenced by the presence 
of transverse pressure-balanced magnetic structures or transverse velocity shears 
[Ghosh et al., 1998]. In this respect it is important to note that in our case the 
mean field direction is constantly changing due to the time-varying contribution of 
magnetic field components to the total magnetic field (compare Figures 5,7). Nev- 
ertheless, the small scale power of the Bz (vertical) fluctuations is the strongest 
all the time. The vertical GSM direction, however, is not always the perpendicular 
direction to the mean magnetic field. Therefore the effect of the mean field, even if 
it represents a robust way for the development of scale-dependent anisotropy has 
to be complemented by other mechanisms, which can produce strong fluctuations 
in perpendicular directions. For example the above mentioned transverse mag- 
netic structures or velocity shears can also redirect the initially parallel-propagating 
Alfven waves to highly oblique waves [Ghosh et al, 1998]. In the plasma sheet, 
BBF-associated dipolarization, increase of Bz (or By when the current sheet is 
tilted), can represent such a transverse magnetic structure. Also, the relative thin- 
ness of the plasma sheet in the z direction and the associated boundary effects can 
produce anisotropy statistics in plasma sheet turbulence [Voros et ah, 2004]. To 
evaluate the possible contributions of these mechanisms a wider statistical study is 
required, however. 

Let us consider now the estimated scaling indices a in comparison with the pre- 
dictions of Goldreich-Sridhar model. When the anisotropy effects are neglected, 
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the conditional statistics explained above gives for BBF-associated flows the large 
scale index a s (js, j'4) = 2.6 ± 0.1, independent on magnetic field component. It is 
definitely different from the Goldreich-Sridhar model which predicts Kolmogorov 
scaling with a ~ 5/3 in perpendicular direction. We could advance an argu- 
ment that the Goldreich-Sridhar anisotropy is weaker at the scales (73, J4), but the 
small-scale a s (j\, J2) reaches the same range of values when the estimated Cj sn 
is large enough and allows to introduce corrections. a s (j3, J4) decreases and ap- 
proaches the Kolmogorov scaling during non-flow or post-BBF intervals, when 
the ratio of the large-scale mean field to the magnetic fluctuations becomes larger 
and the conditions for the development of Goldreich-Sridhar anisotropy are more 
favourable. In this case the estimated a s (js, j'4) = 1.7 ± 0.4. The large error 
appears because the post-flow intervals usually correspond to transient, probably 
decaying turbulence. This also causes that one cannot distinguish between the 
Kolmogorov and Iroshnikov-Kraichnan models. The high value of a s (j3, j'4) dur- 
ing BBF-associated intervals does not necessarily mean that the Goldreich-Sridhar 
phenomenology fails. Magnetic fluctuations are subject to an additional constraint 
which acts during rapid flows only. The limited vertical dimension of the flow 
channel or the plasma sheet thickness can modify the actual value of a increas- 
ing it to a ~ 3, which is a limiting value for fully developed 2D hydrodynamic 
turbulence [Volwerk et al, 2003, 2004, Voros et al, 2004]. 

5 Conclusions 

The work described here demonstrates the ability of the windowed wavelet estima- 
tor for obtaining reliable estimates of the effective second-order scaling parameters 
(cf, a) even in situations when the noise level is relatively high. Magnetic turbu- 
lence seems to be steady only during relatively short time intervals (a few minutes), 
otherwise a exhibits nonstationary behaviour at all scales, which can partially ex- 
plain the large scatter of previously estimated values of scaling indices. As we 
observed magnetic turbulence associated with bursty flows and post-flow intervals, 
our findings agree well with previous knowledge that the flow of plasma in the 
Earth's magnetotail is intrinsically unsteady [e.g. Baumjohann, 2002]. Joint con- 
ditional statistics from two Cluster spacecraft s/c 1 ,3 shows a high level of spatial 
intermittence which is related to the differing scaling features in physically differ- 
ent regions in the plasma sheet. 

The examination of BBF-associated high frequency magnetic fluctuations within 
a sliding window showed that the small-scale (ji, j'2) and large-scale (j'3, j'4) scal- 
ing indices are equal within error estimations. Volwerk et al. [2003, 2004] have 
obtained similar scaling exponents for the time scales 1-13 s. It means that the 
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spectral break predicted by Milovanov et al. [2001] for the scales < 1// C 2 ~ 10 s 
is absent and, using the notation from the Introduction, 012 ~ 03. During non-BBF 
periods the small-scale a.3 ~ a(ji , j'2) aproaches the noise level of the magnetome- 
ter, which introduces an artificial break at the scale ~ Is. Obviously, it cannot be 
interpreted in terms of a new scaling regime describing a scale dependent physics. 

We found that, for the scales analysed, the conversion of the observed scalings 
from frequency to wavenumber representation is model dependent and cannot be 
verified by direct comparison of spatial and temporal variations. It makes the in- 
terpretation of the results in terms of MHD phenomenology more difficult. The 
anisotropy features of magnetic fluctuations provided additional tools to facilitate 
the interpretation of the results. The anisotropy is better visible in the normalized 
power of fluctuations Cf than in the scaling index a. It was shown that large-scale 
magnetic fluctuations are more isotropic while small-scale magnetic turbulence 
evolves in direction perpendicular to the mean magnetic field, usually in the GSM 
Z or Z — Y direction, and the observed anisotropy is scale dependent. It partly 
agrees with the predictions of Goldreich-Sridhar model of MHD turbulence, but 
the additional effect of transverse magnetic structures or velocity shears cannot be 
excluded. On this basis we think the turbulence in the plasma sheet is a mixture of 
Alfvenic wavy turbulence and of eddy 2D turbulence. The latter becomes stronger 
and more intermittent with decreasing scales. 

The occurrence of reconnection and/or rapid tailward flow can increase the 
level of anisotropy. The build-up of a flow channel seems to modify the actual 
value of the scaling indices, indicating the presence of additional constraints on 
Alfvenic and/or eddy turbulence. These are observations which should be verified 
on larger data sets. 

Acknowledgement 

The authors acknowledge the use of the code for the estimation of scaling expo- 
nents developed by P. Abry and his colleagues. 

References 

[1] Abry, P., P. Goncalves, and P. Flandrin, Wavelets, spectrum estimation and 1/f 
processes, Wavelets and Statistics, edited by A. Antoniadis, and G. Oppenheim, 
Lect. Notes Stat., 103, 15, 1995. 

[2] Abry, P., and F. Sellan, Some remarks on and fast implementation for the 
wavelet-based synthesis of fractional Brownian motion proposed by F. Sellan 
and Y. Meyer, Appl. Comput. Harmon. An., 3, 377, 1996. 



20 



[3] Abry, P., P. Flandrin, M.S. Taqqu, and D. Veitch, Wavelets for the analysis, 
estimation and synthesis of scaling data, Self-Similar Network Traffic and Per- 
formance Evaluation, edited by K. Park, and W. Willinger, Wiley Interscience, 
New York, 39, 2000. 

[4] Angelopoulos, V.,W. Baumjohann, C. F. Kennel, F. V. Coroniti, M. G. Kivel- 
son, R. Pellat, R. J. Walker, H. Liihr, and G. Paschmann, Bursty bulk flows in 
the inner plasma sheet, J. Geophys. Res., 97, 4027, 1992. 

[5] Angelopoulos, V., et al., Characteristics of ion flow in the quiet state of the 
inner plasma sheet, Geophys. Res. Lett., 20, 1711, 1993. 

[6] Baker, D.N., et al., Timing of magnetic reconnection during a global 
magnetospheric substorm onset, Geophys. Res. Lett, 29(24), 2190, 
doi:10.1029/2002GL015539, 2002. 

[7] Balogh, A. et al., The Cluster magnetic field investigation: overview of in- 
flight performance and initial results, Ann. Geophys., 19, 1207, 2001. 

[8] Bauer, T. M., W. Baumjohann, R. A. Treumann, N. Sckopke, and H. Liihr, 
Low-frequency waves in the near-Earth plasma sheet, /. Geophys. Res., 100, 
9605, 1995. 

[9] Baumjohann, W., G. Paschmann, and C. A. Cattell, Average plasma properties 
in the central plasma sheet, /. Geophys. Res., 94, 6597, 1989. 

[10] Baumjohann, W., G. Paschmann, and H. Liihr, Characteristics of high-speed 
ion flow in the plasma sheet J. Geophys. Res., 95, 3801, 1990. 

[11] Baumjohann, W., Modes of convection in the magnetotail, Phys. Plasmas, 9, 
3665, 2002. 

[12] Biskamp, D., Magnetohydrodynamic turbulence, Cambridge Univ. Press, 
2003. 

[13] Borovsky, J. E., R. C. Elphic, H. O. Funsten, and M. F. Thomsen, The Earth's 
plasma sheet as a laboratory for flow turbulence in high-/? MHD, /. Plasma 
Phys., 57, 1, 1997. 

[14] Borovsky, J. E., and H. O. Funsten, MHD turbulence in the Earth's plasma 
sheet: Dynamics, dissipation and driving, /. Geophys. Res., 108, 1284, 
doi : 1 . 1 029/2002JA009625 , 2003 . 

[15] Cho, J., and E. T. Vishniac, The anisotropy of magnetohydrodynamic 
Alfvenic turbulence, Astrophys. J., 539, 273, 2000. 



21 



[16] Cho, J., and A., Lazarian, Compressible magnetohydrodynamic turbulence: 
mode coupling, scaling relations, anisotropy, viscosity-damped regime, and as- 
trophysical implications, Mon. Not. R. Astron. Soc, in press, 2004. 

[17] Cho, J., A. Lazarian, and E. T. Vishniac, MHD turbulence: scaling laws and 
astrophysical implications, Lect. Notes Phys., 614, 56, 2003. 

[18] Coroniti, F. V., L. A. Frank, D. J. Williams, R. P. Lepping, F. L. Scarf, S. M. 
Krimigis, and G. Gloeckler, Variability of plasma sheet dynamics, it J. Geophys. 
Res., 85, 2957, 1980. 

[19] Eastman, T. E., L. A. Frank, W. K. Peterson, and W. Lennartsson, The plasma 
sheet boundary layer, /. Geophys. Res., 89, 1553, 1984. 

[20] Goldreich, P., and S. Sridhar, Toward a theory of interstellar turbulence. II. 
Strong Alfvenic turbulence, Astrophys. J., 438, 763, 1995. 

[21] Ghosh, S., W. H. Matthaeus, D. A. Roberts and M. L. Goldstein, The evolu- 
tion of slab fluctuations in the presence of pressure-balanced magnetic structures 
and velocity shears, /. Geophys. Res., 103, 23691, 1998. 

[22] Hoshino, M., A. Nishida, T. Yamamoto, and S. Kokubun, Turbulent mag- 
netic field in the distant magnetotail: Bottom-up process of plasmoid formation? 
Geophys. Res. Lett, 21, 2935, 1994. 

[23] Horbury, T. S., Cluster II analysis of turbulence using correlation functions, 
ESA SP-449, , 89, 2000. 

[24] Katsev, S., and I. L'Heureux, Are Hurst exponents estimated from short or 
irregular time series meaningful?, Comput. Geosci, 29, 1085, 2003. 

[25] Matthaeus, W. H., S. Oughton, S. Ghosh, and M. Hossain, Scaling of 
anisotropy in hydromagnetic turbulence, Phys. Rev. Lett., 81, 2056, 1998. 

[26] Milovanov, A. V., L. M. Zelenyi, G. Zimbardo, and P. Veltri, Self-organized 
branching of magnetotail current systems near the percolation threshold, J. Geo- 
phys. Res., 106, 6291,2001. 

[27] Nakamura, R., et al., Motion of the dipolarization front during a 
flow burst event observed by Cluster, Geophys. Res. Lett, 29(20), 1942, 
doi:10.1029/2002GL015763, 2002. 

[28] Nagai, T, I. Shinohara, M. Fujimoto, S. Machida, R. Nakamura, Y. Saito, and 
T. Mukai, Structure of the Hall current system in the vicinity of the magnetic re- 
connection site,/. Geophys. Res., 108(A10), doi:10.1029/2003JA009900, 2003. 



22 



[29] Neagu, E., J. E. Borovsky, M. F. Thomsen, S. P. Gary, W. Baumjohann, and 
R. A. Treumann, Statistical survey of magnetic field and ion velocity fluctua- 
tions in the near-Earth plasma sheet: Active Magnetospheric Particle Trace Ex- 
plorers/Ion Release Module (AMPTE/IRM) measurements, /. Geophys. Res., 
107(A7), doi:10.1029/2001JA000318, 2002. 

[30] Ohtani, S., T. Higuchi, A. T. Y. Lui, and K. Takahashi, Magnetic fluctuations 
associated with tail current disruption: Fractal analysis, /. Geophys. Res., 100, 
19135, 1995. 

[31] Oughton, S., E. R. Priest and W. H. Matthaeus, The influence of a mean 
magnetic field on three-dimensional magnetohydrodynamic turbulence, /. Fluid 
Mech., 280, 95, 1994. 

[32] Reme, H., et al., First multispacecraft ion measurements in and near the 
Earth's magnetosphere with the identical Cluster ion spectrometry (CIS) ex- 
periment, Ann. Geophys., 19, 1303, 2001. 

[33] Runov, A., et al., Current sheet structure near magnetic X-line observed by 
Cluster, Geophys. Res. Lett, 30(11), 1579, doi:10.1029/2002GL016730, 2003. 

[34] Schodel, R., W. Baumjohann, R. Nakamura, V. A. Seregeev, and T. Mukai, 
Rapid flux transport in the central plasma sheet, /. Geophys. Res., 106, 301, 
2001. 

[35] Sergeev, V. A., V. Angelopoulos, J. T. Gosling, C. A. Cattell, and C. T. Rus- 
sel, Detection of localized, plasma-depleted flux tubes or bubbles in the midtail 
plasma sheet, /. Geophys. Res., 101, 10817, 1996. 

[36] Shebalin, J. V., W. H. Matthaeus, and D. C. Montgomery, Anisotropy in MHD 
turbulence due to a mean magnetic field, J.Plasma Phys., 29, 525, 1983. 

[37] Veitch, D., and P. Abry, A wavelet based joint estimator of the parameters of 
long-range dependence, IEEE Trans. Inf. Theor., 45, 878, 1999. 

[38] Volwerk, M., et al., A statistical study of compressional waves in the tail 
current sheet, Geophys. Res. Lett, 108(A12), 1429, doi:10.1029/2003JA010155, 
2003. 

[39] Volwerk, M., et al., Compressional waves in the neutral sheet, Ann. Geophys., 
22, 303, 2004. 

[40] Voros, Z., et al, Multi-scale magnetic field intermittence in the plasma sheet, 
Ann. Geophys., 21, 1955, 2003. 



23 



[41] Voros, Z., et al., Wavelet analysis of magnetic turbulence in the Earth's 
plasma sheet, Phys. Plasmas, 11, 1333, doi: 10.1063/1.1667499, 2004. 

[42] Watkins, N. W., M. P. Freeman, C. S. Rhodes, and G. Rowlands, Ambiguities 
in determination of self-affinity in the AE-index time series, Fractals, 9, 471, 
2001. 



24 



40 



8 

1 

o 



a. 


MM a 

+ B x 


b. 




a =1.6 

mod 









18:42 18:48 18:54 19:00 19:06 19:12 



Figure 1: Estimation of the scaling index a; (a) Fractal signals A and B to be 
embedded and quiet time Bx measured by Cluster 3 on September 7, 2001; the 
generated fractal signals A and B have the same a mo d = 1.6, but different am- 
plitudes, (b) a a and as from mixed signal; (the dashed line corresponds to the 
theoretical value a mo d = 1.6) 
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Figure 3: The a sn (c'j sn ) relationship (correction curves) for model fractal signals 
with a mo d = 1.5 and 2, computed at the octaves j = and j = 4 
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Figure 4: Estimation of the quiet time noise level of the scaling parameters; a.) Bx 
component of the magnetic field on September 7, 2001 b.) a sn c.) c^ sn ; all for s/c 
1 and 3 
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Figure 5: Two-point (s/c 1,3) estimation of the scaling parameters in comparison 
with the magnetic field and velocity measurements; a.) the Bx components of the 
magnetic field; b.) a sn estimated at small scales (ji, j'2) ~ (0.08, 0.33)s and at 
large scales (73, j'4) = (0.7, 5)s; c.) (cj sn ) estimated at the scale j = 4 ~ 0.33 s; 
the capital letters A and B show two cases when corrections to small-scale a are 
introduced in Figure 6; d.) x-component of perpendicular velocity 
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Figure 6: Correction curves computed for model fractals with a mo d = 
1.5,2,2.5,3; the correction curves reach asymptotically the theoretical values 
o<-mod for Cf sn — > oo and the noise level for Cf sn — > 1; the letter A corresponds 
to s/c 3 while B to s/c 1 ; the error bars at the left-top corner arise from the uncer- 
tainty of parameter estimation in logscale diagram and from the uncertainty of the 
magnetometer noise level estimation 
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Figure 8: Log-log scatterplots of powers at the scale j = 4 ~ 0.33 s; a.)c'j sn (Bx) 
versus Cj sn (5y); b.) c^ sn {Bx) versus c'j sn (Bz) - the crosses correspond to an 
interval of magnetic reconnection 
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Figure 9: Log-log scatterplots of powers at the scale j = 6 ~ 1.3 s; a.)c'j sn (Bx) 
versus Cj sn (i?Y); b.) Cj- sn (Bx) versus c'j sn (Bz) - the crosses correspond to an 
interval of magnetic reconnection 
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Figure 10: Small-scale (ji, j'2) scatterplots of the scaling indices; a.) a sn (Bx) 
versus a sn (B Y ); b.)a sn (B x ) versus a sn (B z ) 
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Figure 11: Magnetic fluctuations in the mean magnetic field coordinate system; 
a.) the perpendicular component (B±); b.) the parallel component (-B11); c.) the 
normalized powers of B± and Bu at the time scale 0.33s; d.) the normalized powers 
of B± and Bit at the time scale 1.3s; 
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Figure 12: Joint conditional statistics in a coordinate system parallel and perpen- 
dicular to the local mean magnetic field at scales (j'3, j'4); a.) scatterplot of all avail- 
able points from s/c 1,3; b.) scatterplot of BBF-associated data points - the crosses 
correspond to an interval of magnetic reconnection; c.) scatterplot of non-BBF- 
associated data points - the crosses correspond to the intervals when \Bx\ > 20 
nT, while the points represents intervals with \Bx\ < 15 nT 
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